Aggregation of superparamagnetic colloids in magnetic fields: the quest for the equilibrium state. 
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Experimental and simulation studies of superparamagnetic colloids in strong external fields have systemati- 
cally shown an irreversible aggregation process in which chains of particles steadily grow and the average size 
increases with time as a power-law. Here we show, by employing Langevin dynamics simulations the existence 
of a different aggregation behavior: aggregates form during a transient period and the system attains an equilib- 
rium distribution of aggregate sizes. A thermodynamic self-assembly theory supports the simulation results and 
it also predicts that the average aggregate size in the equilibrium state depends only on a dimensionless parame- 
ter combining the volume fraction of colloids cf>o and the magnetic coupling parameter V. The conditions under 
which this new behavior can be observed are discussed. 

PACS numbers: 83.10.Mj, 61.43.Hv, 82.70.Dd, 83.80.Gv 



Colloidal aggregation is a subject of active research for both 
practical (e.g. stability of many industrial products) and fun- 
damental reasons (as a test field for statistical-mechanical the- 
ories, for example). Our interest here is in the new physics 
arising in the aggregation behavior of superparamagnetic col- 
loids. These systems are a successful example of implemen- 
tation of a new behavior typical of the nanoscale (superpara- 
magnetism) in new materials with many exciting practical ap- 
plications, ranging from environmental waste capture HI to 
biomedicine Q. Superparamagnetic materials show a large 
magnetic dipole in presence of external field, saturation mag- 
netization similar to that of ferromagnetic materials but no co- 
ercitivity nor remanence at the working temperature. Super- 
paramagnetic colloids are typically made by embedding su- 
perparamagnetic nanocrystals in a non-magnetic matrix (such 
as polystyrene, nanoporous silica or others) [3 1. 

Extensive experimental studies 0^8] as well as computer 
simulations |9] [TO) show that, after application of strong ho- 
mogeneous and inhomogeneous magnetic fields, superparam- 
agnetic colloids form linear aggregates. These chain-like ag- 
gregates increase in length with time, typically with a kinetic 
law compatible with a scale-free, power-law behavior J4][5][8]. 
An important property, typical of dispersions of superparam- 
agnetic colloids, is the reversibility of chain formation: after 
removal of the external magnetic field, the chains rapidly dis- 
aggregate and the initial dispersion (no aggregation) is recov- 
ered 13 El El- Theoretical analysis of experimental results 
have focused on Smoluchowsky rate equations and the appro- 
priate kernels which reproduce the observed kinetics of chain 
growth under applied field J4j [5] . 

In this work, we propose the existence of a different sce- 
nario for superparamagnetic colloids under strong external 
fields. Let us first note that the aggregation of superparam- 
agnetic colloids under external field is, apparently, similar to 
other self-assembly processes such as micelle or gel forma- 
tion. This similarity suggests that chain growth could lead, 
under appropriate conditions, to an equilibrium state with a 
constant mean chain size (and a definite distribution of chains 
of different sizes), as it happens in these other self-assembly 



processes. To the best of our knowledge, this hypothetical 
equilibrium state has never been reported in experiments or 
in simulations, suggesting that it could be difficult to realize 
under the conditions probed in previous studies. In this let- 
ter we will present thermodynamical arguments and brownian 
dynamics simulations supporting the existence of this equi- 
librium state under certain realistic combinations of size and 
saturation magnetization of the colloids. At this point, we 
should stress that understanding the aggregation process of 
superparamagnetic colloids is not only relevant from a fun- 
damental perspective, but it has also practical importance. A 
paradigmatic example is the fast magnetophoretic separation 
process employed in biotechnological applications fUHHEL 
which requires the formation of chains of superparamagnetic 
colloids. 

As in previous simulations (9] [10) we would like to con- 
sider here the minimal model describing superparamagnetic 
colloids: spheres of diameter d with a magnetic dipole dif- 
fusing in a fluid with viscosity rj. For the sake of simplicity, 
we assume here (as in [10 1) that the magnetization of the col- 
loids has reached saturation. This means that each colloid has 
a constant dipole m s (corresponding to saturation magnetiza- 
tion) parallel to the external applied field. This situation is also 
commonly found in experiments (typically at applied fields 
>0.1 T, see ReflTTl [TOtiTIl ). In this situation, the magnetic 
effects can be described by the magnetic coupling parameter 
r defined as the ratio between the maximum of the magnetic 
dipole-dipole attraction and the thermal energy: 
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Hence, our model is characterized by two dimensionless pa- 
rameters, the coupling constant V and the volume fraction of 
colloids (fio. 

The three-dimensional simulations reported here are based 
on a numerical integration of the Langevin stochastic equa- 
tion of motion for each colloid, as in previous works (9] [TO] . 
In this framework, the force acting on each particle is given 
by the sum of a particle-particle interaction force, the viscous 
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FIG. 1: Time evolution of the average size of chains in simulations 
with 0.5 g/1 concentration (4>o = 5.23 x 10 -4 ) and different magnetic 
coupling constants T = 3, 10, 11, 12, 15 and 40. The dashed line is 
power-law fit n ~ t z to the case V = 40. 



FIG. 2: Time evolution of the average chain size obtained in simu- 
lations with r = 10 and four different concentrations (c=0.5, 1, 2.5 
and 5 g/1). Dashed lines indicate the mean aggregate size value at 
equilibrium. 



drag acting against each colloid and a stochastic force cor- 
responding to the thermal noise. The particle-particle inter- 
action potential is given by the sum of the magnetic dipole- 
dipole interaction and an steric, short range strong repulsion 
which prevents overlap between particles. We have neglected 
the effect of sedimentation, considering that our colloids have 
a density of 1 g/cm 3 (which is similar to that of many com- 
mercial superparamagnetic particles since it helps to avoid 
storage problems). All simulations were performed using the 
Langevin dynamics option as implemented in the 21May2008 
version of the LAMMPS program [14|. The equation of mo- 
tion was solved using a time step of 1 ns. Also, we employed 
a very large cutoff (lOd) for the magnetic interactions in order 
to ensure accuracy of the results, although the resulting sim- 
ulations were extremely time consuming and difficult to par- 
allelize. Each second of simulation time requires (depending 
on the specific simulation) between 300-1100 hours of com- 
puter time employing 16 Itanium Montvale processors. Fur- 
ther technical details and movies illustrating the simulations 
are available in the accompanying EPAPS material! 15 1 . In all 
simulations, we had No = 8000 colloids in the simulation 
box, and different concentrations were obtained by employ- 
ing different system volumes V. The diameter of the colloids 
was fixed to d = 100 nm, a value typical for small super- 
paramagetic colloids (although this value is not essential in 
the sense that simulations with the same value of T are ex- 
pected to give equivalent results). We have considered two 
different simulation sets. In the first set (Fig. [T| we consid- 
ered a volume fraction of colloids </>o = 5.23 x 10~ 4 and 
different values of T. In the second simulation set (Fig. |2]i 
we have considered T = 10 and volume fractions cf>o — 
5.23 x 10~ 4 ,1.05 x 10- 3 ,2.62 x 10" 3 and 5.23 x 10~ 3 (which 
corresponds to concentrations between ^0.5 g/L and 5 g/L, 
typical of experiments). 

As a check for the validity of our simulations, we looked 
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FIG. 3: Fraction of aggregates of size s, n s /'^2 s n s in the equilib- 
rium state from simulations with V — 10 and four different con- 
centrations (c=0.5, 1, 2.5 and 5 g/1). Inset: Rescaled distribution of 
aggregates, N a <f>o as a function of s/n. 



for the typical power-law kinetic behavior observed in previ- 
ous works, which consider large values of V (for example in 
the simulations of Ref.OU T is between 100 and 3000). We 
have found (see Fig. [T} that from V as small as T = 15 (and 
4>o = 5.23 x 10~ 4 ) the mean number of colloids per aggregate 
follows h ~ t z with kinetic exponent z « 0.63. Simulations 
with T = 40 and cj) — 5.23 x 10~ 4 give a slightly larger 
kinetic exponent, z ~ 0.645. Our findings are consistent with 
previous works 0|5][TO), which typically report kinetic expo- 
nents in the range 0.6 - 0.7. 

Now we move towards our objective of finding an equilib- 
rium state, with a value of n independent of time. During 
the time scales of our simulations, the equilibrium state was 
reached for simulations with V < 11, as shown in Figures [T] 
and [2] After a transient process of chain growth, we reach a 
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time independent value of h. Note that the size of the aggre- 
gates in the equilibrium state depends strongly on T and 0o- 
At volume fraction <j!>o = 5.23 x 10~ 4 of colloids, we found 
n ~ 3.2 for T = 11 and h ~ 1.7 for T = 10 (see Fig. [T). As 
concentration increases, the equilibrium value for n increases 
from ~ 1.7 at 0.5 g/L to ~ 6.8 at 5 g/L (see Fig. [2). 

Figure [3] displays the distribution of aggregate sizes in the 
equilibrium state (where N s is the number of aggregates of 
size s and n s — N s /V). The main figure shows that the frac- 
tion of aggregates of size s, n s j J2 S n s decays exponentially 
for large s. The inset shows that, after an appropriate nor- 
malization, the fraction of aggregates of size s for different 
concentrations approximately colapses in a single curve. 

The obtained simulation results can be understood by con- 
sidering a simple thermodynamic calculation based on the 
self-assembly theory 1 13 1 originally developed to describe the 
formation of micelles by amphiphilic molecules. In fact, the 
only ingredient which we need to modify in the theory is the 
driving force for micellization (the hydrophobic effect) which 
will be replaced by the magnetic interaction. Let us start by 
considering that the magnetic energy of a chain or aggregate 
made of s dipoles is given by ps — (s — l)e m , i-e. the en- 
ergy arises from s — 1 bonds each one with magnetic energy 
— e m . Although this approximation may seem rather crude, 
the resulting formalism captures the main features of our sim- 
ulation results, as we will see. The chemical potential fi s of a 
colloid which forms part of a chain of s colloids is given by 
the ideal (entropic) term plus the interaction energy term (see 
Eqs.(16.1) and (16.6) in Ref.flj)): 
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where <fi s is the volume fraction of aggregates containing s 
colloids, related to the number of aggregates per unit volume 
through n s ~ <fi s /s. In the equilibrium state, the chemical po- 
tential fii of a colloid in dispersion in a non-aggregated state 
is equal to the chemical potential of a colloid in any of the 
possible chains of size s, so we have /ii — fi s . Using Eq.Q 
we obtain (/3 = l/k B T): 
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By defining x = 0ie^ £m , the constraint 0o = Y^T=i <^s SU P' 
plies, provided x < 1: 

(1 — x) z 

Hence, the x parameter governing the distribution of aggre- 
gate sizes in Eq. ([3]) is controlled by a single quantity defined 

as: 



N* = v / 0oe /3e 



(5) 



In the case N* << 1, Eq. (|4]) yields cj>-y ps (j)Q, i.e. no ag- 
gregation, as expected when the magnetic attraction is weak 
or the system is very diluted. In the opposite case, N* >> 1, 
one finds x ps (1 — l/N*), so that Eq.([3]) supplies: 

n s ~ (1-1/N*) S ^e~ s / N \ (6) 
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FIG. 4: Average number of colloids per aggregate as a function of a 
single parameter combining the volume fraction cj>o and the magnetic 
coupling r. The symbols correspond to simulation results for differ- 
ent combinations of <f>o and T and the dashed line corresponds to the 
theoretical approximation discussed in the text. 



In this approximation, the average size of aggregates gives 
n Ri N* . Also, note that an exponential decay is observed 
in our simulation results for large s (see Fig[3]), thereby sup- 
porting the plausibility of the simplifications introduced in the 
thermodynamic calculation. 

In order to make more explicit predictions, we need to 
relate e m (and hence N*) with known magnetic properties 
of the colloids. To this end, consider an aggregate of two 
colloids in contact, each one with dipole m s (correspond- 
ing to saturation magnetization) parallel to the external ap- 
plied field. Their magnetic interaction energy is given by 
U m = -fc s rr(l-(3/2) sin 2 9) where is the angle between 
the magnetic field and the line joining the centers of the two 
colloids. This interaction is attractive for \0\ < 6q — 54.7°, 
being maximum at 8 = (f3U m (0) = —T). The ther- 
mal average of this interaction over all orientations corre- 
sponding to the bonding between two colloids is given by 
p < U m {6) >= -r[l - | < sin 2 9 >} where: 
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Hence, we have (3 < U m (9) >~ 1 — T. The last equal- 
ity in Eq.|7]) comes from the stationary phase approximation, 
which in this case is extremely good. For example, for T = 10 
we obtain f3 < U m >= —8.96 ~ —9 from a numerical 
evaluation of Eq.(j7]l. Even for smaller values of V the ap- 
proximation is quite good, for example for T = 3 we have 
f3 < U m >= —2.071 fa —2. Therefore, as an estimation for 
e m we take this thermal averaged dipole-dipole interaction, 
which gives f3e m ps T — 1 and then we obtain: 

n^N* ^ v^oe^- 1 ). (8) 
The comparison shown in Fig|4]demonstrates that Eq.([8]l pro- 
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vides a fairly good approximation to the actual average size 
of aggregates observed in simulations. Although the num- 
ber of simulation data points is small and more statistics 
should be desirable, it is remarkable that our simulation re- 
sults are consistent with a universal behavior of the form 
ft = f((f>oe^ r ~ 1 ' 1 ), as expected from the simple thermody- 
namic calculation. The good performance of the model is 
quite remarkable, given its simplicity. The result given by Eq. 
(|8]l could be very useful in practical situations since it can be 
easily evaluated from characterization data (particle size, sat- 
uration magnetization and concentration) measurable in real 
colloidal dispersions. However, in applying Eq.|8]l in a real 
situation, we should keep in mind that our theoretical analy- 
sis is valid only provided that N* is much smaller than the 
number of colloids in the sample, N* << N . 

In view of these analytical results, it is interesting to discuss 
again our results for the simulations which are not observed 
to reach an equilibrium state, at least during our simulation 
times (see Figure [TJ. In the case of V = 40 and <fio — 5.23 x 
1CT 4 , Eq.(|5]) gives N* ~ 6.7 x 10 6 and an estimation from 
Fig{T]suggests an equilibration time of the order of ~10 years. 
Clearly, in this case the reason for the no observation of the 
equilibrium state is that it is unphysical and it will never be 
observed in a real experiment. This argument also justifies 
why the equilibrium state is not observed in previous works, 
since in applications one is typically interested in large values 
of r in order to obtain strong magnetic effects (for example, 
in our previous work we employed V ~ 10 3 J7)). 

The other situation showing a power law in Fig{T](r = 15 
and <f>o — 5.23 x 10 -4 ) corresponds to a very different case. 
Eq.pl) gives iV*=25 and the extrapolation of the kinetics of 
Fig. [TJ suggests an equilibration time of the order of 10 s. 
This calculation suggest that, in this case, we do not observe 
the equilibrium state due to the limitations in computational 
time of the simulations. However, in a real experimental sit- 
uation it should be possible to observe, in this case, an initial 
kinetics obeying a power law followed by the more slower ap- 
proach to an equilibrium state which will contain aggregates 
of substantial sizes. Values for T around 10-15 can be easily 
obtained experimentally by using superparamagnetic colloids 



with d = 100 nm and saturation magnetization of ^30 emu/g 
(see for example 01). Hence, the behavior reported here is 
readily accessible in real lab situations. 
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